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ABSTRACT 

The  temperature  field  produced  by  a  cryosurgical  probe 
embedded  in  a  clear  gelatin-water  test  medium  has  been 
experimentally  investigated.   The  temperature  field  emanating 
from  the  cryoprobe  was  measured  using  a  thin  mylar  sheet 
sprayed  with  liquid  crystals,  a  temperature  sensitive  materi- 
al.  Probe  tip  temperatures   varying  from  -36 °C  to  -117 °C 
were  studied.   The  experimental  results  compared  within 
approximately  9%  of  a  one  dimensional  analytical  solution 
for  predicting  ice  growth  rate. 
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I.   INTRODUCTION 

The  use  of  low  temperatures  as  a  medical  tool  is  not 
new.   As  early  as  3500  B.C.  the  Egyptians  used  ice  as  a  local 
anesthetic  by  applying  it  to  wounds  and  other  ailments. 
However,  it  was  not  until  1851  that  extreme  cold  was  used  as 
a  surgical  procedure.   Dr.  James  Arnott,  a  British  physician, 
used  a  small  container  filled  with  an  ice-brine  mixture  at 
-24°C  as  a  surgical  instrument.   The  container  was  placed  in 
direct  contact  with  skin  cancers  and  was  found  helpful  in  re- 
stricting the  growth  of  the  cancerous  tissue.   In  1961,  with 
the  development  of  the  first  cryosurgical  probe  for  the  treat- 
ment of  Parkinson's  Disease,  cryosurgery  became  an  establish- 
ed surgical  procedure.   Dr.  Irving  Cooper  of  the  Saint  Barna- 
bas Hospital  in  New  York  is  acknowledged  as  the  founder  of 
the  present  day  cryosurgical  technique.   Cooper  developed  a 
small  multi-tubed  probe  with  a  hemispherical  tip.   The  tem- 
perature of  the  tip  could  be  controlled  from  body  temperature, 
37°C,to  -196°C  by  circulating  liquid  nitrogen  through  it 

[41. 

The  cryoprobe  has  primary  and  secondary  applications  in 
treating  diseases  of  the  brain.   The  major  use  has  been  in 
the  production  of  discrete  destructive  thalamic  lesions  in 
the  surgical  treatment  of  Parkinson's  disease  and  related 
disorders.   Unfortunately,  the  application  of  the  cryosur- 
gical technique  to  the  management  of  benign  or  malignant 
tumors  has  been  severely  limited  due  to  an  inability  to 
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predict  lesion  growth  rate  and  ultimate  lesion  size  [17] . 
There  have  been  several  analytical  models  developed  for 
predicting  lesion  growth  rate  and  size  [1,  6,  7].   However, 
these  models  have  yet  to  be  verified  experimentally. 

The  objectives  of  this  thesis  were:   1)   to  experimen- 
tally determine  the  growth  rate  of  the  frozen  region  around 
a  cryoprobe   2)   to  investigate  the  feasibility  of  using 
liquid  crystals,  a  material  that  indicates  temperature 
through  a  change  in  color,  as  a  temperature  sensing  device 
in  the  unfrozen  region   3)   to  compare  the  experimentally 
determined  temperatures  and  ice  growth  rates  with  existing 
analytical  predictions. 

Two  cryoprobes,  0.5  inches  and  0.187  inches  in  diameter, 
were  designed  and  built  along  the  lines  of  the  standard 
Linde  cryoprobe.   The  larger  probe,  which  was  intended  as  a 
prototype,  produced  ice  regions  that  were  approximately 
spherical  in  nature  and  easily  studied.   The  smaller  probe, 
due  to  faulty  design,  produced  irregular  regions  that  did  not 
submit  to  detailed  study.   A  0.15%  gelatin  99.85%  water  test 
medium  was  used  to  simulate  tissue.   A  range  of  probe  tip 
temperatures  varying  from  -36°C  to  -117°C  were  used  to 
create  the  frozen  regions.   In  all  cases  studied,  experimen- 
tal and  analytical  predictions  of  ice  growth  rates  compared 
within  the  estimated  experimental  uncertainty  of  9%.   The 
details  of  the  magnitude  and  extent  of  the  temperature  field 
in  the  unfrozen  region  were  successfully  obtained  with  the 
liquid  crystal  material. 
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II.   DESIGN  OF  THE  SURGICAL  CRYOPROBE 

The  general  surgical  cryoprobe  used  for  necrotizing 
brain  tumors  consists  of  three  concentric  tubes  [17]  .   The 
coolant,  which  is  normally  liquid  nitrogen  (LN2) /  flows 
down  the  inner  tube  and  absorbs  heat  from  the  tip,  causing 
the  coolant  to  vaporize.   The  vapor  flows  back  through  the 
middle  annulus.   The  outer  annulus  is  evacuated  to  minimize 
heat  transfer  between  the  gelatin  and  the  returning  vapor 
(Figure  1) .   Probe  tip  temperature  is  controlled  by  the 
coolant  flow  rate,  which  is  in  turn  controlled  by  regulating 
the  outlet  pressure  from  the  coolant  supply  reservoir. 

The  original  desire  was  to  determine  the  ice  formation 
rate  on  an  existing  Linde  cryosurgical  unit.   Attempts  to 
secure  such  a  unit  were  unsuccessful  resulting  in  a  need  to 
construct  a  unit.   Using  the  basic  design  of  a  Linde  cryo- 
probe [17]  as  a  model,  two  probes  were  designed  and  fabricated 
from  stainless  steel  hypodermic  tubing.   The  first  cryoprobe 
built  was  a  prototype.   It  was  0.5  inches  in  diameter  and 
5.75  inches  in  length  (Figure  2,  3).   The  prototype  was 
tested  and  found  to  work  satisfactorily.   The  same  design 
was  applied  to  a  smaller  cryoprobe,  0.187  inches  in  diameter. 
This  cryoprobe  is  the  same  size  as  the  Linde  Cryoprobe  PR-2 
used  for  brain  tumors.   When  the  smaller  cryoprobe  was  test- 
ed, it  failed  to  form  a  spherical  ice  ball.   instead,  an 
elliptical  ice  ball  formed  indicating  two  problem  areas. 
The  vacuum  volume  was  too  small  and  the  tip  was  not  a 
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hemisphere  but  was  shaped  more  like  a  cylinder.   A  new 
design  should  employ  a  smaller  inner  tube  for  the  LN2  to 
flow  through  and  the  tip  should  be  a  hemisphere  as  in  the 
case  of  the  large  probe.   The  remainder  of  this  thesis  will 
be  concerned  with  the  0.5  inch  cryoprobe. 
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0.5   x    5.75    INCH   CRYOPROBE 


FIGURE    2 
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0.187  x  6.56  INCH  CRYOPROBE 


FIGURE  3 
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III.   EXPERIMENTAL  APPARATUS 

The  experimental  set-up  used  in  the  present  studies 
was  similar  to  the  set-up  used  by  Groff  in  his  studies  of 
the  resistively  heated  surgical  probe  [14].   Briefly,  the 
set-up  consisted  of  a  8  inch  cubical  test  cell,  with  the 
cryoprobe  horizontally  inserted  through  a  small  hole  in 
the  center  of  one  side  of  the  test  cell.   The  liquid  crystal 
material,  which  was  sprayed  on  a  mylar  sheet,  was  placed  in 
the  horizontal  midplane  of  the  cryoprobe.   The  probe  had 
thermocouples  welded  on  the  tip  and  on  the  midpoint  of  the 
stem.   The  test  cell  was  filled  with  a  clear  0.15%  gelatin 
/99.8  5%  water  test  medium,  which  had  the  following  thermal 
properties 

a.  kf  =  .0050  cal/sec/cm/°C 

b.  k   =  .0014   cal/sec/cm/°C 

c.  L   =  79.71    cal/g 

The  cryoprobe  was  connected  to  a  LN2  source.   In  order  to 
record  the  visual  data  presented  by  the  liquid  crystals, 
a  Graflex  Camera  was  placed  directly  above  the  cryoprobe 
tip  (Figure  4) . 

A.   LIQUID  NITROGEN  SYSTEM 

The  LN2  was  stored  in  a  Union  Carbide  CE-8,  which  is  a 
dewar  that  uses  the  pressure  of  the  vaporized  LN2  as  the 
driving  potential  to  force  the  LN2  through  a  regulating 
valve.   The  unit  is  also  equipped  with  a  heater  which  can 
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be  used  to  vaporize  the  LN2  to  increase  the  pressure  as 
necessary.   The  maximum  working  pressure  of  the  system  is 
10.5  psig.   The  regulating  valve  is  a  needle  valve  that 
allows  the  operator  to  regulate  the  flow  of  the  LN2  (Figure 
5). 

To  transfer  the  LN2  from  the  dewar  to  the  cryoprobe  a 
connecting  tube  was  designed  and  built  (Figure  6) .   The  tube 
is  essentially  two  concentric  tubes.   The  LN2  flows  through 
the  inner  annulus  and  the  outer  annulus  is  evacuated.   The 
entire  tube  is  encased  in  neoprene  insulation. 

B.   LIQUID  CRYSTALS 

There  are  three  types  of  liquid  crystals;  smectic, 
nematic  and  cholesteric  [12,  13,  15,  16] .   This  experimental 
study  used  only  the  cholesteric  type.   A  cholesteric  liquid 
crystal  responds  to  changes  in  temperature  by  sequentially 
passing  through  the  complete  visual  spectrum  (red  through 
violet)  in  fractions  or  multi-degrees,  depending  on  which 
cholesterol  esters  comprised  the  formulation.   This  color 
phenomenon  is  reversible.   A  very  important  point  is  that 
at  a  certain  temperature  a  given  material  or  combination  of 
materials  will  always  exhibit  the  same  color.   The  light 
scattered  by  the  liquid  crystals  represents  only  a  fraction 
of  the  incident  light.   The  remaining  portion  of  the  incident 
light  is  transmitted  by  the  liquid  crystals.   Therefore,  an 
absorptive  black  background  must  be  used  to  prevent  reflec- 
tion of  the  transmitted  light,  thereby  enhancing  the  color 
resolution. 
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Since  the  colors  exhibited  by  the  cholesteric  liquid 
crystal  are  unique  for  a  specific  temperature,  the  quanti- 
tative measurement  of  temperature  is  possible  to  an  accuracy 
of  approximately  0.1°C/  while  the  speed  of  response  of  the 
liquid  crystals  is  less  than  0.2  seconds. 

C.   DESCRIPTION  OF  LIQUID  CRYSTAL  SHEETS 

Mylar  sheeting  was  chosen  for  the  foundation  of  the 
liquid  crystal,  because  it  is  thin,  pliable  and  has  a  small 
thermal  conductivity,  which  for  this  system  results  in 
negligible  heat  transfer  along  the  sheet. 

The  mylar  sheet  was  cut  into  a  7  15/16"  square  with  a 
4  1/8" x  1/2"  notch  cut  into  the  sheet  for  the  experimental 
runs  (Figure  7) .   The  liquid  crystals  were  applied  to  the 
mylar  sheets  with  a  Paasche,  Model  H  "3  in  1"  air  brush. 
The  steps  taken  in  preparing  the  liquid  crystal  sheets  were: 

A.  Thoroughly  clean  both  sides  of  mylar  sheet  with 
acetone. 

B.  Use  liquid  slurry  received  from  NCR  without  dilution, 

C.  Spray  first  coat  of  liquid  crystals  on  one  side  of 
the  clean  mylar  surface  allowing  the  crystal  to 
dry  for  4  0  to  60  minutes. 

D.  Spray  second  coat  and  allow  to  dry  for  20  minutes. 

E.  Spray  third  coat  and  allow  to  dry  for  one  hour. 

F.  Apply  the  first  coat  of  Testor's  flat  black  paint 
on  the  dryed  liquid  crystals  and  allow  to  dry  for 
30  minutes. 
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LIQUID  CRYSTALS  SPFAYED  ON  A  MYLAP 

SHEET,  THE  SHEET  IS  ENCASED  IN  THE 

LIQUID  CRYSTAL  FRAME. 


FIGURE  7 


22 


G.   Apply  the  second  coat  and  let  dry  for  one  hour. 

H.   Apply  two  coats  of  Rex  Polyurethane  Interior- 
Exterior  77-5  over  the  black  paint  paying  particular 
attention  to  the  edges.   This  will  water  proof  the 
liquid  crystal  sheet. 

The  mylar  sheet  for  the  experimental  run  was  divided 
into  five  areas  (Figure  8).   The  W-5,  W-10,  R-15  and  R-19 
areas  were  covered  with  masking  tape.   W-0  crystals  were 
applied  as  described  earlier,  but  only  one  coat  of  Testor's 
flat  black  paint  was  used.   The  W-5  area  was  uncovered  next 
and  the  W-5  crystals  were  sprayed.   Then  one  coat  Of  Testor'; 
flat  black  paint  was  applied  to  the  W-5  area.   The  process 
was  repeated  for  the  W-10,  R-15  and  R-19  areas.   A  second 
coat  of  Testor's  flat  black  paint  v/as  then  applied  to  the 
entire  sheet  and  allowed  to  dry.   The  sheet  was  then  water 
proofed.   The  reader  is  referred  to  Appendix  A  for  calibra- 
tion technique  and  data. 

D.   TEST  CELL 

The  test  cell  was  composed  of  two  parts.   The  first  was 
the  liquid  crystal  sheet  frame.   The  frame,  which  supported 
the  liquid  crystal  sheet,  was  7  15/16  inches  square  and  was 
made  up  of  two  plexiglass  squares  that  were  screwed  together 
with  the  liquid  crystal  sheet  in  the  middle  (Figure  8) .  On 
one  side  of  the  frame  was  a  0.5  inch  diameter  hole  through 
which  the  cryoprobe  fit. 

The  second  part  was  a  clear  plexiglass  8  inch  cube. 
Four  inches  down  from  the  top  of  the  cube  was  a  0.5  inch 
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wide  shelf  that  supported  the  liquid  crystal  frame  (Figure  9) 
A  0.5  inch  diameter  hole  in  the  side  of  the  box  corresponded 
to  the  hole  in  the  frame.   When  the  cryoprobe  was  placed 
through  the  hole  in  the  box  and  the  frame,  the  longitudinal 
axis  of  the  cryoprobe  was  in  the  same  plane  as  the  liquid 
crystal  sheet. 
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EIGHT  INCH  CUBICAL  TEST  CELL 


FIGURE  9 
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IV.   EXPERIMENTAL  PROCEDURE 

The  liquid  crystal  frame  was  placed  in  the  box.   The 
cryoprobe  was  pushed  through  the  hole  in  the  box  and  the 
frame.   A  copper/constantan  thermocouple  was  welded  on  the 
tip  of  the  cryoprobe.   To  record  the  temperatures,  a  Honey- 
well Electronik  194  dual  strip  recorder  with  ice  junction 
was  used. 

The  gelatin  was  made  and  allowed  to  cool  for  six  hours 
(Appendix  C) .   After  six  hours  the  gelatin  was  cool  enough 
to  pour  into  the  assembled  box  without  thermally  shocking 
the  liquid  crystals.   The  gelatin  then  was  allowed  to  come 
to  room  temperature.   The  LN2  dewar  was  connected  to  the 
cryoprobe  by  means  of  an  insulated  double  wall  tube  (Figure 
6) .   A  vacuum  pump  was  connected  to  the  cryoprobe  and  the 
double  wall  tubing.   A  pressure  of  approximately  10mm  of 
mercury  was  maintained  in  both  probe  and  tubing.   The  Graf lex 
camera  (Appendix  B)  was  positioned  directly  over  the  cryo- 
probe and  all  necessary  adjustments  were  made.   The  ther- 
mocouples were  connected  to  the  strip  recorder.   The  pressure 
on  the  LN2  dewar  was  between  3  and  9.5  psi. 

At  time  equal  zero  the  LN2  was  turned  on.   The  tip 
temperature  was  monitored  with  the  recorder.   When  the  tip 
temperature  reached  a  predetermined  temperature,  the  controls 
on  the  dewar  were  adjusted  so  the  tip  temperature  remained 
constant.   Pictures  were  taken  at  times  of  0,  1,  3,  5,  10, 
15,  30  and  every  15  minutes  thereafter  until  the  21.5°C 
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isotherm  moved  to  the  edge  of  the  box.   At  this  point  in 
time  the  experiment  was  terminated. 

The  photographs  recorded  the  isotherms  and  ice  ball 
diameter.   Using  photo  interpolation  (knowing  the  probe 
diameter  to  be  0.50  inches)  the  actual  ice  ball  diameter 
and  the  distance  from  the  cryoprobe  tip  center  to  the 
various  isotherms  were  determined. 

Six  experimental  runs  were  made  varying  the  probe  tip 
temperatures  as  shown  below: 

Experiment        Tip  Temperature       Initial  Gelatin 

Temperature 

1.  -36°C  22°C 

2.  -53°C  21°C 

3.  -65°C  21°C 

4.  -75°C  20°C 

5.  -97°C  25.5°C 

6.  -117°C  22.5°C 
In  all  cases  the  phase  change  temperature  was  -2.0°C. 

This  was  found  by  allowing  the  strip  chart  recorder  to  run 
for  one  hour  after  the  experiment  ended.   The  probe  tip 
thermocouple  steadied  on  -2°C  and  remained  there  until  the 
ice  ball  had  completely  melted  away  from  the  cryoprobe. 

Figure  10  shows  a  segment  of  chart  paper  for  a  typical 
experimental  run.   The  first  half  is  the  response  when  the 
LN2  is  turned  on  and  the  second  half  is  the  response  when 
the  LN2  is  turned  off. 

Figures  11,  12,  13,  and  14  show  the  ice  ball  forming  in 
the  gelatin  with  the  liquid  crystals  visually  presenting 
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the  isotherms. 

The 

rel 

ationship 

between  the 

colors  and 

temperatures  are: 

Red 

Green 

Blue 

W-0 

1.7°C 

6.2°C 

8.0°C 

W-5 

3.3°C 

9.3°C 

11.0°C 

W-10 

10.0°C 

12.7°C 

15.4°C 

R-15 

15.9°C 

16.5°C 

17.1°C 

R-19 

19.8°C 

20.8°C 

21.5°C 
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0.5  INCH  CRYOPROBE  EMBEDDED  IN  THE 
GELATIN  SURROUNDED  BY  THE  LIQUID  CRYSTAL 
SHEET  AT  TIME  EQUAL  ZERO.   SEE  FIGURE  7 
FOR  COLOR  TEMPERATURE  RELATIONSHIPS. 


FIGURE  11 
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0.5  INCH  CRYOPROBE  EMBEDDED  IN  THE 
GELATIN  SURROUNDED  BY  THE  LIQUID 
CRYSTAL  SHEET  AT  TIME  EQUAL  FIVE  MINUTES. 
SEE  FIGURE  7  FOR  COLOR  TEMPERATURE  RELATIONSHIPS 


FIGURE  12 
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0.5  INCH  CRYOPROBE  EMBEDDED  IN  THE  GELATIN      • 
SURROUNDED  BY  THE  LIQUID  CRYSTAL 
SHEET  AT  TIME  EQUAL  THIRTY  MINUTES 
SEE  FIGURE  7  FOR  COLOR  TEMPERATURE  RELATIONSHIPS, 


FIGURE  13 
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0.5  INCH  CRYOPROBE  EMBEDDED  IN  THE  GELATIN 

SURROUNDED  BY  THE  LIQUID  CRYSTAL  SHEET 
AT  TIME  EQUAL  SIXTY  MINUTES.   SEE  FIGURE  7 
FOR  COLOR  TEMPERATURE  RELATIONSHIPS. 


FIGURE  14 
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V.   THEORETICAL  ANALYSIS 

When  a  cryoprobe  is  inserted  into  tissue  and  brought  to 
a  cryogenic  temperature,  a  thermal  field  develops  in  the 
tissue.   A  lesion  is  formed  as  the  tissue  in  the  immediate 
vicinity  of  the  cryoprobe  undergoes  a  phase  change.   The 
freezing  front  will  propagate  through  the  tissue  until  steady 
state  is  reached  within  the  tissue.   The  specific  heat,  den- 
sity, thermal  conductivity  and  latent  heat  of  fusion,  as  well 
as  blood  flow  and  metabolism  all  contribute  in  controlling 
ultimate  lesion  size  and  growth  rate  for  a  given  cryoprobe 
diameter  and  particular  tip  temperature.   Therefore,  to 
predict  lesion  size  in  a  structure  such  as  a  human  brain, 
it  is  necessary  to  know  not  only  the  cryoprobe  tip  temperature, 
but  also  the  thermal  properties  of  the  structure  immediately 
adjacent  to  the  cryoprobe  tip  [4] . 

There  have  been  several  models  [1,  6,  7]  that  describe 
the  transient  behavior  of  a  cryolesion.   Cooper  and  Trezek 
[6,  7]  present  an  approximate  solution  to  the  problem  of 
transient  freezing  around  cylinders  and  spheres  embedded  in 
tissue.   Their  model  included  terms  for  the  thermal  proper- 
ties, as  well  as  for  blood  flow  and  metabolism.   While  the 
blood  flow  and  metabolism  terms  are  important,  they  have 
been  disregarded  in  this  study.   In  order  to  take  the 
phenomenon  into  account,  a  live  brain  is  necessary.   Since 
the  school  lacked  the  proper  facilities  for  operating  on 
experimental  animals,  it  was  decided  to  experimentally 
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investigate  the  rate  of  ice  formation  produced  by  a  cryo- 
probe  embedded  in  gelatin.   The  experimental  results  were 
compared  to  the  analytical  solution  of  Cooper  and  Trezek 
to  ascertain  if  their  basic  premise  (neglecting  heat  capacity 
effects  in  both  the  frozen  and  unfrozen  phases)  is  valid. 

The  following  is  a  brief  development  of  the  solution 
presented  by  Cooper  and  Trezek  without  metabolism  and  blood 
flow  terms. 

The  energy  equation  (see  Nomenclature  for  definition  of 
terms) 

PC   II  =  kV2T  (1) 

is  used  to  describe  the  temperature  field  around  the  cryo- 
probe.   The  energy  equation  will  take  the  following  form 
when  used  to  describe  one  dimensional  heat  transfer  in  the 
frozen  and  unfrozen  regions  surrounding  a  spherical  probe. 
Unfrozen  region 


pc 

3T 

at 

:  kV2T  =  i  J_   (r2 
r2   3r   \ 

3T  ] 
3r  / 

(2) 

Frozen 

region 

3T 
pCf  at 

=  kfV2Tf  =  H*   J-   ( 

f    f    r2   3r   \ 

3T 
r2  — 1 
3r, 

) 

(3) 

l  the  following  boundary  conditions 

Tf  =  Ts 

(4) 

r  = 

Tf  =  Tpc  =  T 

(5) 

r  -*-  * 

T   =  T 

o 

(6) 

r  = 

k    Hi  =    k  n 

Kf   3r       K   3r 

+  pL 

3K 

3t 

(7) 
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By  introducing  the  following  non-dimensional  groups 


e  -Tf~ 

To 

T 

o 

T   -  T 

ft  -         ° 
6    T   -  T 
s     o 

T    -  T 
0    _  pc    o 

pc    T    -  T 
r     s     o 

«~t 

r*  =  2- 

o 

_  k(TQ  -  Tpc)t 

PLro2 

the  normalized  energy  equation  becomes: 
frozen  region 

R2   9R    \    3R  /        kf   L  ^ 


unfrozen  region 

A.   _L   fR2  ii  \  =  °  (To  "  TPC}    il  (9) 

R*    3r   \   9R/         L        3T 

with  the  boundary  conditions  4,  5,  6,  and  7  becoming 

R  =  1  ( 

R  =  r*  f 

R   ->   00  ( 

36f 
R  =  r*     kf  W"   =k  3R 

3  T 

Analytical  solutions  of  equation  (8)  and  (9)  satisfying 

conditions  (10),  (11),  (12),  and  (13)  are  not  known.  However, 

k  cf  (T   -  T   )     c  (T   -  T   ) 
if  we  assume  r-^ 2£  and  ^= £°—  are  small,  that 

is,  if  the  heat  capacity  effects  are  small  compared  to  latent 

heat  effects  (8)  and  (9)  take  the  approximate  form: 

frozen  region 

1   3  L,  36f 


f =  1 

(10) 

f  =  0pc  = 

=  e 

(11) 

->  0 

(12) 

14-  ke 

.  3r* 

(13) 

,R2  — -  )=  0  (14) 

3R    \    3R 


37 


unfrozen  region 


19  l=  0  (15) 


R*     3R   \     9R 
the  solutions  to  equations  (14)  and  (15)  are 
unfrozen  region 


R 
frozen  region 


(16) 


(£)( 


Vj  (1  ~  > 


Equations  (18)  and  (19)  are  substituted  into  equation  (13) 
to  find  N 

ar*  AVi  V l  -y   )    i  (18) 

dT   \k  /f^/v  +  (r*"i)/    r* 

Further  rearrangement  of  equation  (18)  yields 

dx  =  r*  (1  "  r*)    dr*  (19) 

(r*+n  ) 

where 


kf 

pc 


dr* 
at  steady  state  ^ —  =  0,   this  necessitates  from  equation 

(19)  that 

r*  =  -  n  (20) 

at  steady  state.   Integrating  (19)  between  the  limits  1  and 

-n  yields: 

x  =  (n+D  (r*-l)  +  n(n+l)ln^i  -  \   (r*-l)   (21) 

Using  equation  (20)  and  (21)  an  analytical  prediction 
of  the  rate  of  growth  of  the  ice  ball  can  be  made.   A 
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FORTRAN  program  (Appendix  E)  was  written  that  uses  equation 
(20)  and  (21)  to  predict  the  rate  of  lesion  growth.   The 
program  also  predicts  time  to  steady  state  and  maximum 
lesion  size  for  the  particular  cryoprobe  size  and  tip  tem- 
perature. 

When  using  the  FORTRAN  program  written  in  conjunction 
with  this  thesis  the  actual  radius  of  the  hemispherical 
cryoprobe  tip  is  used  vice  a  modified  radius.   At  the 
beginning  of  this  study  it  was  felt  that  perhaps  the  radius 
should  be  the  radius  of  an  imaginary  sphere  that  has  the 
same  surface  area  as  the  hemispherical  cryoprobe  tip  since 
the  theoretical  analysis  assumed  a   spherical  case.   However 
the  best  agreement  between  the  theoretical  solution  and 
experimental  data  occurred  when  the  actual  hemispherical 
radius  of  the  cryoprobe  tip  was  used  in  the  FORTRAN  program. 
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VI.   RESULTS 

Correlation  between  the  experimental  data  and  the  theory 
was  excellent.   The  theoretical  curves  for  the  six  different 
cryoprobe  tip  temperatures,  as  calculated  by  the  FORTRAN 
program  (Appendix  E) ,  are  illustrated  in  Figures  15,  16,  17, 
18,  19,  and  20.   The  experimental  data  is  superimposed  on 
the  graphs  to  demonstrate  the  agreement.   For  the  first 
hour  the  correlation  is  almost  exact,  then  the  experimental 
data  starts  to  rise  above  the  theoretical  curve,  but  remains 
within  the  uncertainty  bounds  (Appendix  D) .   The  rise  in  the 
experimental  data  can  be  explained  by  the  theoretical 
boundary  conditions  not  being  met  after  the  first  hour. 
That  is,  the  19°C  isotherm  had  moved  to  the  edge  of  the  test 
cell  and  the  gelatin  no  longer  qualified  as  an  "infinite" 
medium. 

Figure  21  is  a  typical  graph  of  the  nondimensional 
time  vs.  nondimensional  radius  for  a  tip  temperature  of 
-25°C,   The  solid  line  represents  the  theoretical  solution 
and  the  dashed  lines  represent  the  uncertainty  bounds. 
The  reader  is  referred  to  Appendix  D  for  details  on  the  uncer- 
tainty analysis. 

The  transient  temperature  field  shows  the  details  of 
a  lesion  formation.   An  example  is  illustrated  in  Figure 
22.   The  discontinuity  at  -2°C,  the  phase  change  temperature, 
indicates  the  lesion  size  at  a  particular  time  after  initiat- 
ing cooling.   The  experimental  points  on  the  graph  were 
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COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
PREDICTION  OF  ICE  FRONT  LOCATION,  r*,  AS  A 
FUNCTION  OF  TIME,  t,  FOR  A  PROBE  TEMPERATURE 
OF  -36°C. 

FIGURE  15 
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COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  PRE- 
DICTION OF  ICE  FRONT  LOCATION,  r*,  AS  A  FUNCTION 
OF  TIME,  t,  FOR  A  PROBE  TEMPERATURE  OF  -53 °C. 


FIGURE  16 
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o   Experimentally  determined  data 
—  Theory 
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COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  PRE- 
DICTION OF  ICE  FRONT  LOCATION,  r*,  AS  A  FUNCTION 
OF  TIME, t,  FOR  A  PROBE  TEMPERATURE  OF  -65°C. 


FIGURE  17 
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COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  PRE- 
DICTION OF  ICE  FRONT  LOCATION,  r*,  AS  A  FUNCTION 
OF  TIME,  x,  FOR  A  PROBE  TEMPERATURE  OF  -75°C. 


FIGURE  18 
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COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  PRE- 
DICTION OF  ICE  FRONT  LOCATION,  r* ,    AS  A  FUNCTION 
OF  TIME,  t,  FOR  A  PROBE  TEMPERATURE  OF  -97 °C. 


FIGURE  19 
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COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  PRE- 
DICTION OF  ICE  FRONT  LOCATION,  r*,  AS  A  FUNCTION 
OF  TIME,  x,  FOR  A  PROBE  TEMPERATURE  OF  -177 °C. 


FIGURE  20 
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GRAPH  OF  THE  UNCERTAINTY  BOUNDS  OF  THE  THEORET- 
ICAL CURVE  AND  UNCERTAINTY  LIMITS  OF  THE  EXPERI- 
MENTALLY DETERMINED  DATA  AT  A  CRYOPROBE  TIP 
TEMPERATURE.  OF  -7  5°C. 

FIGURE  21 
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o   Liquid  Crystal  Data,  Time  =  5  mins . 

o   Liquid  Crystal  Data,  Time  =  30  mins. 

x   Liquid  Crystal  Data,  Time  =  90  mins. 

Liquid  Crystal  Data,  Time  =  210  mins 
—  Theory 
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GRAPH  OF  THE  TRANSIENT  TEMPERATURE  FIELD  AT 
TIMES  OF  5,  30,  90,  210  MINUTES  AT  A  CRYOPROBE 
TIP  TEMPERATURE  OF  -75°C. 

FIGURE  22 
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determined  from  the  liquid  crystal  sheet.   The  experimentally 
determined  points  show  that  the  temperature  gradient  in  the 
unfrozen  phase  is  much  sharper  than  the  theoretically  pre- 
dicted temperature  gradient.   This  difference  is  attributed 
to  the  heat  capacity  effect  which  is  neglected  in  the  theo- 
retical analysis.   However,  the  thermal  capacitance  does 
not  seem  to  be  important  in  determining  the  ice  growth  rate. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  use  of  liquid  crystals  afforded  an  unique  opportu- 
nity to  visually  witness  the  rate  of  movement  of  various 
isotherms.   The  analytical  prediction  of  the  lesion  growth 
rate,  as  presented  by  Cooper  and  Trezek,  shows  good  agree- 
ment with  experimentally  determined  data.   The  assumption 
made  in  the  theoretical  analysis,  that  the  heat  capacity 
in  the  frozen  phase  is  negligible,  seems  reasonable.   How- 
ever, the  thermal  capacity  in  the  unfrozen  region  does  have 
a  definite  influence  on  the  extent  of  the  temperature  field. 

The  ability  to  control  the  tip  temperature  was  difficult 
at  higher  temperatures  (-36°C,  -53°C)  and  relatively  easy 
at  lower  temperatures  (-117 °C) .   It  is  recommended  that  an 
automatic  control  system  be  developed  to  control  the  cryo- 
probe  tip  temperature. 

The  smaller  cryoprobe  should  function  properly  if  the 
tip  were  redesigned  so  as  to  be  hemispherical.   Also  when  the 
small  cryoprobe  is  used  in  an  experimental  run,  low  LN2 
pressure  (1.5  to  3.5  psig)  should  be  used. 
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APPENDIX  A 

Liquid  Crystals 
For  this  project  seven  different  temperature  ranges  of 
the  cholesteric  type  were  used.   The  manufacturer  was 
National  Cash  Register  (NCR) .   The  crystals  used  were: 
(NCR  Notation) 

A.  W-0 

B.  W-5 

C.  W-10 

D.  R-15 

E.  R-19 

F.  R-23 

G.  R-27 

W  -  wide  temperature  range  -5°C 

R  -  regular  temperature  range  4°C 

No  information  was  received  as  to  the  exact  temperature  at 

which  color  changes  would  occur  nor  the  exact  meaning  of 

the  ranges  5°/4°C.   Thus  it  was  necessary  to  calibrate  the 

crystals  prior  to  any  experimentation  to  find  their  ranges 

of  utilization. 

Two  calibration  mylar  sheets  were  prepared.   Both  sheets 
were  four  inch  square.   The  same  technique  used  for  preparing 
the  liquid  crystal  sheet  used  in  the  experiment  was  used  for 
preparing  the  calibration  sheets.   Each  sheet  was  divided 
into  fourths.   The  first  sheet  contained  the  W-0,  W-5, 
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W-10,  and  R-15  crystals.   The  second  sheet  was  similarly 
divided  with  R-19,  R-23,  and  R-27  applied.   Each  sheet  was 
calibrated  separately  in  the  Rosemount,  model  913A,  calibra- 
tion bath  (constant  temperature  bath) .   The  bath  was  lowered 
to  0°C.   The  first  sheet,  containing  W-0,  W-5,  W-10  and 
R-15,  was  placed  in  the  bath.   The  bath  temperature  was 
monitored  by  a  Rosemount  Commutating  Bridge,  model  92A,  and 
galvanometer  arrangement  with  a  Rosemount  Platinum  Resistance 
Temperature  Standard  162C.   Throughout  the  calibration  pro- 
cedure one  observer  read  the  color  changes  while  the  other 
read  the  commutating  Bridge.   The  following  results  were 
obtained. 

Calibration  of  Liquid  Crystals 


W-0 

W-5 

W-10 

R-15 

R-19 

R-23 

R-27 


Red 

Green 

Blue 

1.7°C 

6.2°C 

8.0°C 

3.3°C 

9.3°C 

11.0°C 

10.0°C 

12.7°C 

15.4°C 

15.9°C 

16.5°C 

17.1°C 

19.8°C 

20.8°C 

21.5°C 

22.4°C 

23.5°C 

24.6°C 

26.7°C 

27.9°C 

28.9°C 
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APPENDIX  B 

Photographic  Techniques 

In  order  to  make  a  permanent  record  of  the  experiments, 
pictures  were  required.   For  data  taking  Polaroid  Polacolor 
4x5  land  film  Packets,  color/ type  58,  were  used  in  a  Graf lex 
camera  (4x5  press  Graflex) .   The  experimental  runs  were 
made  at  night.   The  lighting  was  supplied  by  a  Tungsten 
3200°K  light.   The  film  used  was  a  daylight  type  necessitat- 
ing the  use  of  a  Tiffen  Photar  Filter  80B  series  6  for  color. 
When  taking  pictures  for  publication  quality,  35mm  high 
speed  Ektachrome  color  film  was  used  in  a  Honeywell  Pentax 
H-3.   This  film  was  designed  for  use  with  the  Tungsten 
3200°K  lights. 

In  both  cases  the  technique  used  for  taking  the  pictures 
was  the  same.   The  camera,  on  a  tripod,  was  placed  directly 
over  the  cryoprobe.   The  camera  was  focused  and  settings 
were  made.   The  Tungsten  light  was  placed  approximately  six 
feet  away  from  the  cryoprobe.   When  a  picture  was  desired, 
the  Tungsten  3200°K  light  was  turned  on,  the  room  lights 
turned  off  and  the  picture  taken.   As  soon  as  the  picture 
was  taken,  the  Tungsten  light  was  turned  off.   The  total 
time  the  Tungsten  light  was  left  on  was  approximately  ten 
seconds.   It  was  felt  that  the  heat  generated  by  the  light 
had  little  or  no  effect  on  the  freezing  rate. 

The  general  photographic  techniques  were  used  when 
taking  pictures  through  the  gelatin.   The  highest  f  stop 
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with  the  slowest  shutter  speed  was  used  to  get  the  greatest 
depth  of  field.   The  main  drawback  with  the  Polaroid  film 
was  that  it  was  a  daylight  film  with  a  reasonable  ASA,  but 
when  the  filter  was  added  the  ASA  was  lowered  to  24.   Color 
reproduction  was  considered  adequate  for  data  taking. 
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APPENDIX  C 

Preparing  the  Gelatin 
To  make  the  gelatin  the  following  recipe  was  used: 

a.  Heat  6800  ml  of  water  to  boiling 

b.  Soak  120  gms  of  Schilling's  Gelatin  in  1200 

ml  of  cold  water  for  3  minutes.   Do  not  disturb 
during  this  time. 

c.  Pour  boiling  water  into  the  gelatin  mixture. 

d.  Stir  gently  while  pouring  in  the  hot  water. 

If  stirred  too  vigorously  a  foam  will  appear  on 
the  surface. 

e.  After  all  hot  water  is  in  gelatin  mixture,  stir 
to  insure  all  gelatin  is  thoroughly  mixed. 

f .  Leave  gelatin  in  mixing  container  and  put 

in  refrigerator  (initial  temperature  32°F)  for 
3  hours. 

g.  Take  out  of  refrigerator  and  let  stand  at  room 
temperature  for  3  hours. 

h.   Pour  into  experimental  box. 

i.   Let  stand  at  room  temperature  for  12  hours. 
Gelatin  is  now  ready  for  use. 
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APPENDIX  D 


The  uncertainty  in  a  function  f  (x,y,z)  computed  from 
measured  values  of  several  independent  quantities  (x,y,z) 
may  be  estimated  as 


^vpr^T*  (if  ¥f  '*  (sTJ 


(1) 


where  wf  is  the  estimated  uncertainty  in  f  due  to  estimated 

uncertainties  of  u  ,  w  ,  w   in  x,  y,  z  respectively  [2,  11] 

Such  an  uncertainty  analysis  was  carried  out  on  the 
k   /T    -  T  \  k  (T   -  T   )  t 

n  =  TT  t^ T'     ~1  and  T  =  2 E2 

7  pLro2 

to  ascertain  the  influence  of  experimental  uncertainties  on 
the  theoretical  and  experimental  determinations  of  r* . 
Applying  (1)  to  n  and  x  it  is  found: 


€^)  •frMW-fe*) 


pc     s  '     \  pc 


and 


The  following  are  the  estimated  true  values  and  uncer- 
tainties for  the  experimental  run  with  the  probe  temperature 
held  at  -75°C. 

T    -  T   =  -77°C  ±  1°C  20:1 

pc    s 

T    -  T   =  -22°C  ±  1°C  20:1 

pc    o 
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k     =      0.0014    ±    .0001    ~i  20:1 

C   cm   sec 


kf    =    0.0050    ±    .0001   £^-  20:1 

C   cm   sec 


p      =   1.00    ±    .01   2H  20:1 

Cm3 

L      =   79.71    ±    .01   ^i-  20:1 

gm 

rQ   =    0.025    ±    .002   Cm  20:1 


WT        -   T 

T  PC.  m   S      -   0.01 
pc  s 


WT         -    T 


&■ -2.      =    0.05 

-    T 
pc  o 


X     =    °'07 


t|-  °-02 


p   _ 
p 

0. 

.01 

L 

0. 

01 

wr 
o 

=    0.01 

The  results  of  the  uncertainty  analysis  indicated  that 
for  this  particular  run: 


— -  =0.09      and     -2.  =0.09 
t  n 

Similar  calculations  were  made  on  the  other  five  experi- 

0J  CO 

mental  runs  and  in  all  cases  —  and  —  fell  between  0.08  and 

t     n 

0.09.   A  typical  uncertainty  plot  is  shown  in  Figure  21. 
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